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DLFFERENTIAL VARIATIONS IN BALLISTICS, WITH 
APPLICATIONS TO THE QUALITATIVE PROPER- 
TIES OF THE TRAJECTORY * 

By 
T. H. GRONWALL 

The investigation of the influence on the trajectory of changes in initial 
velocity, atmospheric conditions, etc., is most conveniently performed by the 
method of differential variations. These satisfy a system of linear differential 
equations of the fourth order, which was first set up in its most general form 
by F. R. Moulton. By the systematic use of the system of differential 
equations adjoint to the preceding one. Bliss has developed a very simple and 
practical method for computing all the differential variations.! 

This method makes use of the existence of a first integral X = const, (for 
explanation of notations, see section 1) by means of which the adjoint system 
reduces from the fourth to the third order. In the present paper, another 
first integral x' X + y' n + x" v + y" p = const, is established, and the ad- 
joint system thereby reduced to the second order (sections 1-3), and in 
consequence, the numerical computation of the variations is materially 
shortened. 

Of greater theoretical interest is the fact that, the system of linear differential 
equations involved being of the second order, the general behavior of their 
solutions may be determined in a fairly complete manner. This being done, 
a large number of qualitative properties of the trajectory follow almost imme- 
diately (sections 4 and 5). In the last section, attention is called to some 
unsolved problems of the same kind. 

1. The diffeeential equations for the variations, and the adjoint 

SYSTEM 

The retardation (in the direction of the tangent to the trajectory) of a pro- 
jectile moving with velocity v through air of normal density is taken to be 
vG{v)/C, where the positive function G (v) depends on v alone and is the 

* Presented to the Society, Oct. 25, 1919. 

t G. A. Bliss: I. A method of computing differential corrections for a trajectory, Journal 
U. S. Artillery, vol. 51 (1919), p. 445. II. Differential equations containing arbitrary 
functions, these Transactions, vol. 21 (1920), pp. 79-92. III. Functions of lines in 
ballistics, these Transactions, vol. 21 (1920), pp. 93-106. 
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same for all projectiles, while the ballistic coefBcient C depends on the weight 
and shape of the projectile. The decrease in atmospheric density with 
increasing altitude y is taken into account by multiplying the previous expres- 
sion for the retardation by a positive factor H (y) where H (0) = 1 .* 

Referring the projectile to a right-handed system of rectangular coordinates 
X, y ,z with the origin at the muzzle of the gun, the positive a;-axis horizontal 
and directed toward the target, the positive y-axis vertical and directed 
upward, so that the positive z-axis is directed to the right of the line of fire, 
the differential equations of the trajectory aref 

x" = -Ex', 

(1) y" = -Ey' - g, 

z" = - Ez', 

where g is the acceleration of gravity, 

(2) j,^Gi.)Hiy) 

and 

(3) V = VZ+TN^'" 
The initial conditions are for i = , 

a;o = , x'o = Vo cos a , 

(4) 2/0 = 0, 2/0 = ■Wo sin a, 

So = 0, Zo = 0, 

where «o is the initial velocity and a the angle of departure. Since 
{z'/x'Y = by (1), it follows from (4) that z = 0. 

Let us now consider the variations in this trajectory introduced by changes 
bvo in the initial velocity and ha in the angle of departure, or by changes in the 
retardation due to wind, to changes in the functions G{v) and H {y) produced 
by other atmospheric conditions, and to a change in the ballistic coefficient C . 
DenotebyGi(«) = G{v) + 8G{v),Hi{y) = H (y) + SH(y),Ci = C + 8C 
the new G- and H-functions and the new ballistic coefficient, and the com- 
ponents of the wind velocity w, which we assume to depend on the altitude y 
alone, hy w^ = w^iy) , Wy = Wy{y) , w^ = Wz{y) . Let the coordinates of 
the projectile in the changed trajectory he x, y ,z &t the time t; the retarda- 
tion of the projectile, at a given altitude y , depends only on its velocity V 

* In practice, H (y) = e'^-"^^'^^ (lengths being measured in meters and time in seconds), 
while G(v) is given in tabular form. 

t The notations used, which differ in some respects from those of Bliss, are those officially 
adopted by the Technical Staff. 
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relative to the atmosphere, the components of which are x' — w^iy) , 
y' — Wy{y) , z' — w^iy) , so that the expression for the retardation now 
becomes 

(5) vMI)ihiy} = VE.. 

Ci 
Consequently, the differential equations of the changed t ajectory become 
x" = - Erix' -w^y)), 

(6) y" = -E^{y' -wAy))-g, 

z" = - £,{-z' - w^y)) , 
with the initial conditions 

Xa = 0, Xq = {% + 6«o) cos (a + ha) , 

(7) yo = 0, y'o = {vq + hv«) sin (a + ha), 

2o = , 2o = . 
Now write 

(8) l^x = X — X, Ay = y — y , Az = S — z 

and subtract equations (1) from (6); we obtain 

(Ax)" = -E[iAxy -wAy + ^y)] 

- (El - E)[x' + (Ax)' - wAy + Ay)], 

(9) (Ay)" = -E[{Ayy -wAy + Ay)] 

- (E, - E)[y' + (Ay)' - Wy{y + Ay)], 

(Az)"= -E[{Azy -w,{y + Ay)] 

-iE^-E)[z'+iAzy -wAy + Ay)], 

with the initial conditions for < = 

(Aa;)o = 0, (Aa;)d = (vo + 5?)o) cos (a + da) — ■»o cos a, 

(10) iAy)o = 0, (A2/)o = (% + hvo) sin (a + Sa) - Vo sin a, 

(As)o = 0, (A2);=0. 

We shall assume that, for all values of <, we have 

\w^(y + Ay) - w^iy)\ <cW^-Ay, 

where c is a constant and Wx equals the maximum value of Wx{y) along the 
trajectory,* with similar conditions on Wy and w^, and moreover that G {v) 

* This assumption is introduced in order to make the difference w^ {y + Ay) — w^iy) 
a small quantity of the second order in the expansions below. 
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and H(y) have bounded derivatives of the second order. We now expand 
the differences to the right in (9) by Taylor's theorem with remainder term 
of the second order, and denote by dots any terms containing as factors 
squares or products of the following quantities: Ax, Ay, Az, (Ax)' , (Ay)' , 
{Azy,wAy),iVy{y),wAy),8G{v),8G'{v)[ = G[M - G' (v)], dH (y) , 
8H' (y) and SC. We thus obtain 

(11) (Ax)' -w^iy + Ay) = (Ax)' - w^y) + ••• 

and similar expressions in y and z; moreover, writing Wx, Wy, Wz for w^iy) , 
i^y{y) J ^2 ( 2/ ) and observing that z' = , 

F2 = [x' + {Ax)' -wAy + ^y)f 

+ [y' + {Ay)' -Wy{y + Ay)]^ 

+ {{Az)' - w,{y + Ay )]^ 
= i? + 2x'{{Ax)' -w,]+2y'[{Ay)' -Wy]+---; 

V^v + '^liAx)' -w.] + y-[{Ay)'-Wy]+---; 



(12) 



Gi(F) = G,{v)+ {V -t)G[{v) + •■■ 

= G {v) + hG {v) + {V - v)G' {v) + ■ ■ ■ ; 

H,{y)=H,{y)+Aym{y)+ ■■• 

= H{y) + 8H{y) + AyH'{y) + •••; 

Ci C + dC C\ C J' 

G^{V)H^{y) _ G(v)H{y) r SGjv) .y_.G'{v) 
Ci C L Giv) ^^ G(v) 

^ H{y) ^ y H{y) C^ 
and replacing V — ?) by (12) 

w,]+y'[iAyy-Wy]G'{v) 



£,-E = £{^_l(^^ 



V G{v) 

^^^^ +AH'{y),5GM 8Hiy) SC 

y H{y) '^ G{v) '^ H{y) C^ 

We now introduce (11) and (12) in (9), replace Ax, Ay and Az by ^, »? and f 
and omit all the terms denoted by dots; writing ^' = |i and 17' = »7i , we thus 
obtain the differential equations 
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1 I l-nG , nf 1 , liG \ /Til^G , SH dC\ 

(15) r= -£(r' -«'.), 

with the initial conditions for i = 0, obtained from (7) by expanding and 
omitting powers of 5% and da higher than the first, 

1(0) = 0, 
^(0)=0, 

(16) li ( ) = cos a-Svo — Vo sin a • So: = a;o — - — y'o da , 

iji ( ) = sin a • 5?)o + ^o cos a- 8a = Vq— + x'o da, 

(17) r(0)=0, r(0)=0. 

The quantities ^, »?, f defined by (14)-(17) are the differential variations 
(of the first order) of a; , y and z . This appellation is justified by the following 
property, for the proof of which see Bliss III : 

Let € be any positive quantity, and suppose that dvo, da, dC , w^, Wy, Wz, 
dO , dG' , dH and dH' are less than e in absolute value (the last seven at every 
point of the trajectory). Then the expressions 

Ax-^ jAxY - li Ay-r, {Ay)' - rji Az - t (Az)' - ^' 

, , , , , 

€€€€€€ 

tend toward zero with e , and this uniformly in any time interval ^ t ^ to. 

In integrating the equations for the differential variations, we begin with 
(15) and (17), which define the lateral deflection ^ (t) of the projectile at the 
time t . 

From the first of (1), we have E = — x"/x' , and (15) becomes 

i^'/x'Y = Ew,/x', 
whence by (17) 

"'Ew, 



Jo ^ 



510 T. H. GRONWALL [October 

integrating again, we obtain 

(18) r(o = CE{T) ^^*^~''^^K .dT, 

Jo ^T 

as is readily verified by differentiation, and in the particular case where 
W2 = f or r < ^0 , but w^ = const. 5^ for r S #0 > we have 

Integrating by parts, we find the 

Deflection due to a constant cross wind Wz from the time to onward 



(19) 



»'> = [' -^-^^^^i^]- 



To obtain expressions for |, 57, §1 and J71, we follow Bliss in introducing the 
system of linear differential equations adjoint to the system (14), namely 

X' = 0, 

^^^^ /= -\ + Ev + x'E^ix'v + y'p), 

vG 

p' = -fi + Ep + y'E^ix'v + y'p). 

By the fundamental property of the adjoint system, any solution ^ , j? , $1 , r]i 
of (14) and any solution X, fj., v, p oi (20) are connected by the relation 

{^1^ + m + ^i^ + m pY = \ E (1 + x"^\w, + x' y'E^Wy 

+ p,'£g^. + E(l + /g)«,„ 
,^(hG , bH hCX'X 

or reducing the right-hand member by means of the last two equations (20) 
(21) {i\ + w + ^iv + nip)' = (X + /)w, + {ix + p')Wy 

v, ' I , .(bG ,m bc\ 
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Integrating from t = to t = to and using (16), we find 

( ^X + 1JJU + ^1 J' + i7ip)(=(„ = {x' V + y' p)(,~+ (x' p — y' v)oda 

Vo 

(22) + r\\ + v')w,dt+ r\iJL + p')Wydt 

Jo Jo 

-j Eix. + yp)^-+---yt. 

This formula gives us the value of any linear combination of ^, ?;, fi and J71 
at the time t = to, provided that we have obtained a solution of (20) taking 
as initial values for t = to the coefficients of the linear combination in question. 

2. First integrals of the adjoint system, and reduction of the latter 
from the fourth to the second order 

The existence of two first integrals of (20) is due to the fact that the equa- 
tions (1) contain neither x nor t explicitly. Since the equations (1) do not 
contain x explicitly, they are also satisfied by x = x + e, y = y , z = z, 
where e is a constant, or in other words, these x, y, S satisfy the particular 
equations (6) in which Wx = Wy = w^ = dG = 8H = 8C = . We have here 
Ax = €, Ay = (Ax)' = {AyY = 0, and the first order terms in respect to e 
in these differences are $ = «, '7 = ?i = '7i = 0. These quantities must 
satisfy equations (14) (specialized of course to w^, = • • • = 6C = ) , and 
consequently (21) speciahzed in the same way gives for any solution X, fi, v , p 
of (20) upon division by e, (1-X + 0-^ + O-i' + 0-p)' = orX' = 0, whence 
the first integral of (20) 

(23) X = const. 

This integral might obviously have been read off at once from the first of (20), 
but the derivation given emphasizes the analogy to the integral which we shall 
now obtain. Since equations (1) do not contain t explicitly, they are also satis- 
fied hy X = X (t + e) , y = y (t + e) , z = z (t + e) with a constant e , and ex- 
panding the differences Aa; = x{t + e) — x(t) , (Ax)' = x' (t + e) — x' (t) , 
etc., by Taylor's theorem, with remainder term of the second order, we see that 
the first order terms in e are ^ = x' (t)e, tj = y' (t) e , ^1 = x" (t) e and 
Tji = y" {t)e. These must satisfy the system (14) specialized as above, 
whence it follows as before from (21) that for every solution X, ix, v, p oi (20), 
we have (x' "K + y' [x + x" v + y" p)' = or 

(24) K + x''K + y' n + x" V + y" p = (k = const.) 

which is another first integral of (20).* This integral may of course be verified 

by differentiating (24), substituting the values of X' , p,' , v' , p' from (20) and 

* It seems plausible that in general, that is, without specializing G ( v ) or H (y) , there 
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using the following formulas, obtained by differentiating (1) and (2), 

,^^, IdE G' , , H' , „vG' G' , , H' , 

the last expression following from m' = x' x" + y' y" = — Ei? — gy' . To 
perform the reduction of (20) to the second order by means of the two integrals 
(23) and (24), we observe that (1) allows us to write the last two equations 
(20) in the forms 

{x'vy=-x'-K + x'^E^{x'v + y' p), 

vG 

iy'v)' = -y'\-gv + x'y'E^ix'f + y'p), 

(26) % 

(x'p) =-x'ix +x'y'E^{x'v + y'p) , 

C 

iy' pY = - y'lJ' - 9P + y' E—{x' V + y' p), 

whence 

{x' V + y' pY = - x'\ - y'y. - gp + E —■ {x' v + y' p) . 

On the other hand, the substitution of x" and y" from (1) in (24) gives 

(27) gp = K + x'-K + y' iJi- E{x'v + y' p). 
Introducing the notation 

(28) h = h{y) = -^^, 

■n {y) 

and writing down the second of (20) together with the result of substituting 
gp from (27) in the preceding expression for (x' i' + ?/' p)', we find 

m' = -hEix'v + y'p), 

(29) (x'v + y'pY = (l+''-^)E(x'v + y'p) - 2{x'\ + y' n) - k, 

which is the required second order system in n and x' v -\- y' p . Having 

exists no further integral of (20) algebraic m.\, n, v, p,x,y and their derivatives. In special 
cases, however, new integrals may be found. Thus, ior H (y) = const, (corresponding to 
an atmosphere of constant density) the second equation (20) gives n = const., while in the 
case G (v) = const, (so that the retardation at normal air density is proportional to the 
velocity) the third equation (20) combined with (23) gives x' v + x\ = const. 
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integrated (29), 


we obtain p from (27) and v from the identity: 


(30) 


v=^{x'v + y' p) -Kp* 

X X 
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Before proceeding further, we note that (26) gives 

(31) {x'p-y'vY = y'\-x' n + gv 

and that, by means of (1) and (25), the second equation (29) may be written 
in the form 

[x'E{x'v + y'p)]' 

(32) = -{h + g^\ x'y'Eix'v + y'p) - x'E{2x'\ + 2y' n + k) . 

3. Formulas for the variations in range and maximum ordinate 

In the following, coordinates, velocities and times pertaining to the point 
of fall (where 2/ = ) will be denoted by the subscript co , and those pertaining 
to the summit (where ?/' = ) by the subscript s . 

On the trajectory defined hj x = x{t) , y = y {t) the time of flight t^ is 
the positive root of the equation 

(33) 2/(C)=0, 
and the range x^ is given by 

(34) x^ = x{tj.^ 

The angle of fall w is the angle between and ir/2 defined by 

(35) tan CO = - ^ = - ^^ • 

The time ^ required to reach the summit of the trajectory is determined by 
the equation 

(36) 2/'(4)=0 

and the maximum ordinate is then 

(37) ys-yih). 

In the changed trajectory defined by^ = x{t) -\- Kx{t) ,y = y (t) + Ay (t) , 

* The system (29) is adapted to the theoretical purposes of the present paper. For the 
application of the method to the actual computation of differential variations, it is preferable 
to integrate numerically the differential equation of the second order in fi obtained by eliminat- 
ing x' V -\- y' p in (29). A detailed exposition of this method will be published elsewhere. 

t For a proof of the existence of a unique positive root of (33), under the sole assumption 
that E > Q everywhere, see T. H. Gronwall, Qualilalive properties of the ballistic trajectory, 
Annals of Mathematics, ser. II, vol. 21 (1920), pp. 44-65. 
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I = Az(t) , the time of flight t^ + At^ is given by 

(38) y{t^ + AtJ+Ayit^ + AtJ=0. 

Now y {to, + Ato, ) = y (to,) + y' (ti) At^, , where ti lies between t^, and 
to, + Ato,, and since y (tj =0, (25) gives At^, = -Ayit^ + AtJ/y'iti). 
Hence Af„ is of the order of magnitude of Ay , that i > of j; , since the difference 
Ay — r] IS, of higher order of magnitude than 77; consequently Ay{t^ + At^,) 

— vi^u) + (a term of higher order of magnitude), and the first order term in 
Ato,, or the differential variation dt^, in time of flight is given by 

(39) 8t^=-^l^.* 

y (to,) 

The change in range Ax„ is given by 

Xo, + Aa;„ = x{to, + A<„) + Ax{to, + A<„) 
= x{tJ+x'{tJAt^ + i{tJ 

+ [ix"{t,){AtJ^ + Ax(t^ + AtJ -Hi.)] 

with ti between to, and <„ + At^, , the bracketed terms being of higher order 
of magnitude than the first. Retaining only terms of the first order of 
magnitude, we find the differential variation SXo, in range to be x' (to,) 8to, 
+ ^ {to,) , or substituting 6<„ from (39) and using (35) 

(40) 5cr„ = n4) +»?(<.) cot CO. 
The change in maximum ordinate Ay^ is given by 

ys + Ay, = y{ts + At, ) + Ay (ts + Ats) , 

and since 2/' ( 4 ) = , the differences y {t, + At,) — y (ts) and Ay{ts + At, ) 

— ri {ts) are of higher order of magnitude than the first, so that we obtain for 
the differential variation in maximum ordinate 

(41) Sys = ri{ts). 

After these preliminaries, let us return to the remark at the end of section 2. 
By (40), 8X0, equals the linear combination ^-1 + ??• cot co + ?i-0 + tji-G at f 
= to, , and to obtain 8X0, , we have therefore to integrate (20) with the initial 
conditions 

(42) X = 1 , ;u = cot w , V = 0, p = for t = to,. 
Substitution of < = /„ in (23) therefore gives X = 1 , and in (24), 

K + x'o, + y'o, cot w = 

* For the general theorem of which this determination of Sto, is a special application, see 
BUss II, pp. 90-92. 
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or K = by (35). Consequently, for the determination of the differential 
variations in range, the system (29) becomes 

m' = -hE{x'v + y'p), 

(43) {x'v+ y'pY = (^1 +^-^^E{x' V + y' p) - 2(x' + y' ,x) 
with the initial conditions 

(44) ;u = cot CO , x' V + y' p = for t = t^. 
Moreover, (27) becomes 

(45) gp = x' + y'^-E{x'v + y'p), 
and (22), making to = t^ and using (40) 

5x^ = (x' V + y' p)o— + (x' p — y' v)o 8a 

(46) + f" {l + p')w^dt+ f {fj. + p')wydt 

Jo Jo 

Considering separately the various terms to the right in (46), we find the 

following expressions for: 

Range variation due to change 8vo in initial velocity: 



(47) 8x^ = ix' V + y' p)o 



8vo 



Range variation due to change 8a in angle of departure : 
(48) 8x^ = (x' p — y' v)o 8a. 

Range variation due to a following wind Wx'. 



8x^=1 (1 + v')wxdt, 
Jo 



and when w^, = for < < <o, w^ = const. 5^ for < ^ t^: 

Range variation due to a constant following wind Wx from the time 4 onward: 

(49) 8x^ = [t^ — to — V {ta)] Wx 

(since v{t^) = by (42)), and similarly 

Range variation due to a constant vertical wind Wy from the time to onward: 



(50) 8x, 



= 1 iidt - p{to) \wy. 
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Range variation due to a change 8H {y) in H (y): 

E{x'v + y' p)~dt. 
Range variation due to a change 5(7 in the ballistic coefficient: 

(52) ^"^"^^j E^^'^ + y' P)di- 

Passing to the variations in maximum ordinate, (41) shows that hy^ equals 
the linear combination ^-O + ij-l + ^i-O + t/i-O at t = U, and we have to 
integrate (20) with the initial conditions 

(53) X = 0, M = l, v = Q, p = for t = U-* 

Substitution oi t = ts in (23) therefore gives X = 0, and in (24), k = 0, since 
y' (ts) = . Consequently, for the determination of the differential variations 
in maximum ordinate, the system (29) becomes 

m' = ~hE{x'v + y' p), 

(54) ix' V + y' pY = (^l + '^^E{x' V + y' p) -2y' p., 
with the initial conditions 

(55) M = 1 , x' v + y' p = Q for t = t,. 

Moreover, (27) becomes 

(56) gp = y'fx- E{x'p + y' p), 

and (22), making to — ts and using (41) 

Sys = (x' V + y' p)o—+ {x' p - y' v)o 5a 

(57) + f'v'w^dt+ f' {tJi + p')wydt 



dt. 



Separating the terms to the right, we find for the 

Maximum ordinate variation due to change dvo in initial velocity: 

(58) Sys = {x'p + y'p)o—' 



* The two sets of solutions of (20) belonging to range and maximum ordinate variations 
are thus entirely different, since they satisfy different initial conditions (42) and (53). This 
should be constantly borne in mind in the following, as we have refrained from distinguishing 
the two sets by subscripts in order to simplify the typography. 



19213 DIFFERENTIAL VAKIATIONS IN BALLISTICS 517 

Maximum ordinate variation due to change 5a in angle of departure: 

(59) by, = {x' p - y' v)o5a, 

Maximum ordinate variation due to a constant following wind Wx from the 
time to onward: 

(60) 8ys = - v{to)wx, 

Maximum ordinate variation due to a constant vertical wind Wy from the 
time to onward: 

(61) §2/8= I ixdt - p {to) \wy , 
Maximum ordinate variation due to a change SH (y) in H {y): 

E{x'v + y'p)'-^dt, 
Maximum ordinate variation due to a change 5C in the ballistic coefficient: 

(63) ^y^^^J Eix'v + y'p)dt. 

4. Qualitative peopeeties of the differential variations in range 

In addition to assuming G{v) > 0, H (y) >0 as we have done from the 
outset in accordance with the physical significance of these functions, we shall 
introduce, in this and the following sections, the two hypotheses 

(64) G'M>0, hiy)^0, 

the latter being equivalent to H' {y) ^0 by (28) . Physically, this means 
that the retardation at normal density increases more rapidly than the first 
power of the velocity, and that the density of the atmosphere nowhere increases 
with the altitude. 

The fundamental qualitative property of the range variations is expressed 
by the inequality 

(65) x'v + y'p>0 for O^t <t„. 

To prove this, we write 

and obtain from the second equation (43) and (44) 

(66) e*{x'v + y' p) =2 fe'i'ix' + y' ix)dt. 

Jo 

Since {x' + y' p,)' = — E (x' + y' p,) — gp — y' p by (1), and since &it = t„ 
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we have x' + y' n = hy (35), ju = cot co and ;u' = by (43) and (44), it 
follows that {x' + y' ij,y = — g cot w at t = t„ , so that {x' + y' ix)' < for 
t„ — e ^ t ^t„, where e is sufficiently small. Thus x' + y' n decreases 
toward zero in this interval, and (66) shows that x' v + y' p is positive for 
t„ — e < t < t„ , or more generally, in the interior of any time interval ending 
at t = t^ , and in which x' + y' ix ^ 0. Now suppose that x' v + y' p is not 
always positive in the interval ^ t < t^ , then there exists a zero ti nearest 
to t„ so that x' V + y' p = ior t = h but x' v + y' p > ior h < t < t^ , 
and from what precedes it is seen that we must have 

(67) x' + y' IX <0 for t = h. 

The first of (43) shows that ju' ^ for <i < i < /„ , and from (44), we con- 
sequently obtain 

(68) ju S cot CO 

for ti < t < t„ at least. Since a;' > for all values of t* we conclude from 
(67) and (68) that y' (h) < , that is, the hypothetical zero oi x' v + y' p 
must occur, if at all, on the falling branch of the trajectory. For h < t < t„ , 
or more generally, as long as y' < and x' v + y' p > , it follows from the 
third of (26), the fundamental hypothesis (64), and (68) that 

(x' pY < — x' IX ^ — x' cot 03, 

and integrating this inequality from t to <„ and observing that p = for 
t = t„hy (42), we find 

(69) x' p > {x^ — x) cot w 

for ti^t<t„. But from (45) we obtain gp = x' + y' p. for t = h, since 
x' V + y' p = ior this value of t, and in consequence of (67), p must be nega- 
tive att = ti, while (69) shows that p is positive for this value of t , and this 
contradiction proves (65). 

Moreover, it is now seen that (68) is true for ^ < < <„ , the equality sign 
to be taken only when h is identically zero, while (69) holds at least when 
y' ^0, that is, at least for ts = t < t„ . 

From the first of (43) and (65), we have n' ^ ior ^ t < t„ , the equality 
sign to be taken only when h = 0, and consequently, ju decreases toward 
cot CO as < increases from to <„ (except when h vanishes identically, ix being 
then equal to cot co ) . 

We shall now prove, moreover, that x' + y' p. is positive and decreases to zero 
when t increases from to t„. On the rising branch of the trajectory where 
y' > 0, this is evident, since x' and y' are positive and decrease with t increasing, 

* This follows at once from Xj = vo coaa > 0, since (1) gives x' = x„ eJo ■^<" . 
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while fx is positive and decreases or remains constant, li x' -\- y' n < any- 
where on the falling branch, we must have x' + y' ix = and (x' + y' n)' ^ 
at some point on the falling branch, since it was shown that x' + y' ix > 
suflBciently near t = t^. But at such a point — gfx — y' fx' ^ on account of 

{x' +y'iiy = -E{x' +y'ix) ~ gfx - y' fx' , 

and this is impossible since fx > 0, 2/' = 0, m' = 0. Consequently 
x' + y' fi > on the falling branch also (except at the point of fall, where it 
vanishes), and if x' + y' fx does not decrease steadily on the falling branch, 
there must be a minimum, distinct from the point of fall, and at the minimum 
point, {x' + y' fi)' = or — E {x' + y' fi) — gn — y' ix' = , which is im- 
possible, since the first two terms are negative, and the third negative or zero. 

It is now possible to establish a number of qualitative properties of the 
differential variations in range. 

Since, by (65), x' v + y' p > at t = 0, it follows from (47) that 
The range increases when the initial velocity is increased. 

Integrating the second of (43) from t to <„ and using (44), we find 

x'p + y'p= p^2(x' + y'fx) - ^ 1+ ^'^£(a;' . + / p) ]^<, 

or integrating the term in y' ix by parts and using the first of (43), 

x' V + y' p = 2{x^ — x) - 2yn 
(70) „ rf^^^G' _oi,Aj,^^',.^..'.^dt. 



-J"'(l+f-'-2%)^(a;'^ + 2/'p)< 



If we assume that 

(71) 2%<1+^' 

over the entire trajectory, the last term in (70) is negative by (65), and 
making < = so that a; = t/ = 0, we find {x' v + y' p)o < 2a;„ or by (47), 
assuming dvo to be positive, 

(72) ^"<2^», 

X^ Do 

that is 

For low trajectories, where y is so small that (71) is satisfied everywhere, the range 

increases less rapidly than the square of the initial velocity. 

When h = 0, i.e., the atmospheric density is constant, (71) is always true 

since G' > 0; in this particular case, the theorem was proved in an entirely 

different manner by Petrini.* 

* H. Petrini, Om hallisiiska egenskaper hos kastbanor, Arkiv for Matematik och 
F y s i k (Stockholm), vol. 7 (1912), pp. 1-29. 
Trans. Am. Math. Soo. 34. 
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In actual practice, we have vG'/G > 0.51 and h = 0.0001036, so that (71) 
is satisfied whenever the maximum ordinate does not exceed 7500 meters. 
In vacuum, the well-known formula x^ = vl sin 2a/g shows that the range 
increases exactly as the square of the velocity. 

From the first of (26) it follows, since X = 1 by (42) and x' v + y' p > 
by (65), that {x' v + x)' > , so that x' v + x increases with t, and since 
V = Oatt = t^hy (42), we find 

(73) p <^^^ for 0^t<t^. 

x 

Now (31) gives, by means of (68) and (73), since X = 1 and a;' > 0, 

[x p - y v) < y - X cot w + g -, — ; 

integrating from t to i„ and observing (42), we find 

- {x' p - y' v)< J'"l- x' coto, + y' + g'^^^ldt 

r v'T" 

= \ — X cot (jl — (x^ — x)^ \ ) 
L a; J, 

the integration in finite terms being easily verified by means of (1), and finally 

(74) x' p-y' v> {x^-x)( cot o^-'^-^- 

Making i = 0, and comparing to (48), we obtain for the range variation due 
to an increase in the angle of departure 

(75) 5a;„ > a;,, ( cot w — tan a) ha, 

so that bxa > when tan a ^ cot w or a -|- w = ir/2; from the well-known 

fact that a < w , it follows that 

An increase in the angle of departure increases the range when the sum of the 

angles of departure and fall does not exceed 90°, or more particularly, when the 

angle of fall does not exceed 45°. 

This theorem gives no information as to whether the maximum range at a 

given initial velocity occurs at an angle of departure equal to 45° (which is the 

case in vacuum), or smaller or greater. In actual practice, all three cases are 

found to occur, according to the values of the initial velocity and the ballistic 

coefficient.* From (73), we obtain 

t^-to-v{to) > t^ - t, - ^^ ~ ''}^\ 

X {to) 

* A theoretical investigation of the special case of a homogeneous atmosphere ( A = ) 
and a retardation proportional to the rath power of the velocity is given by Petrini, 1. c. 
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and making Wx = Wz in (49) and (19), it follows that 

The increase in range due to a constant following wind from the time to onward 
is greater than the deflection at the point of fall due to a constant cross wind of the 
same velocity from the same time to onward. 
Using (65) in (52), we see that 

An increase in the ballistic coefficient increases the range. 

Similarly it is shown that an increase in atmospheric density ( 5H > every- 
where on the trajectory) or in retardation (5G > 0) will decrease the range. 

5. Qualitative properties of the differential variations in maximum 

ordinate* 

The fundamental inequality for the maximum ordinate variations is 

(76) x'v + y'p>0 for 0Si<4, 

and its proof is simpler than in the case of the range, since now y' does not 
change its sign in the time interval considered. Writing 

we obtain from the second equation (54) and (55) 

(77) e'l' {x' V + y' p) = 2 C ' e* y' ]i.dt 

so that, since y' > on the rising branch of the trajectory, and /* = 1 at 
t = ts, yL\& positive, and consequently also x' v ■\- y' phy (77), for t less than 
but sufficiently close to 4 . If (76) is not true everywhere, there consequently 
exists a ti such that x' v + y' p is zero at ii , but positive for ti < t < tg. 
Then the first of (54) shows that p, decreases or is constant as t increases from 
ti to ts so that ju ^ 1 for <i < i < 4 by (55), and (77) therefore gives a positive 
x' V + y' p at t = ti, contrary to our assumption. Thus (76) is proved. 
From the first of (54) and (76) it follows that ju' ^ so that p. decreases 
steadily toward unity as t increases from to t, (or is constant = 1 when h 
vanishes identically). Using (76), the first of (26) shows, X being zero by (53) 

* The following will explain why no discussion is given of the properties of the variations 
in time of flight, which are more important from a practical point of view than those in maxi- 
mum ordinate. It follows from (39) that the speoiaUzation of (29) proper to the time of flight 
variations is k = 1, X = 0, m=— 1/2/Jj, i' = /> = at t = iu. The non-homogeneous 
terms in the second equation (29) then become — (2y' li + 1) taking the values + 1 at tu 
and — 1 at <,. On account of this change of sign, the method used in the text to prove 
x' V + y' p to be of constant sign is inapplicable. Moreover, numerical computation shows 
the behavior of the time of flight variations to be so much more compUcated and irregular 
than in the cases of range and maximum ordinate, that it appears difficult even to formulate 
conjectures in respect to their quaUtative properties. 
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and (23), that (a;' »')' > 0, so that x' v increases toward the value zero at U, 
and consequently 

(78) f < for ^ < < <s . 

Since y' > 0, the comparison of (76) and (78) shows that 

(79) P>0 for O^Kts, 

and the first of (54), (56) and (79) now give the inequality 

hy' IX + / = gph S 0, 
or by (28) {ix/H (y))' S whence, observing that fj, = 1 atts, 

..^H{y) 



(80) 



'H{y,) 



the equality sign holding only when h is identically zero, or H (y) = const. 
From (54) and (55), we find 

x'v + y'p= £^2y't,-(^l+'^^E{x'v + y'p)yt 

= 2(2/. -2//X)- P'(^l+~-2hy^Eix'p + y'p)dt 

and from (76) and (71) it is seen that for i = we have < {x' v + y' p)o <2ys, 
so that (58) gives for an increase 8vo in initial velocity 

(81) 0<^^<2^ 



or 



For an increase in the initial velocity, the maximum ordinate increases, and for 
low trajectories, where y is so small that (71) is satisfied everywhere, the maximum 
ordinate increases less rapidly than the square of the initial velocity. 
When h = 0, the latter part of the theorem is always true. Comparing (31) 
and (78), it is seen that 

{x' p — y' v)' < — x' p, -^ — x' 

since m = 1 , whence integrating from t to t^ 

(82) x' p - y' V > x^- X 

and making t = and using (59) : 

An increase in the angle of departure increases the maximum ordinate, the increase 

satisfying the inequality 

(83) 8ys>x,8a. 
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The trajectory being concave downward, we evidently have t/s < Xs tan a , 

whence (83) gives 

8ys ^ . J S ( sin a ) 
-^ > cot a 5a = . — - 
ys sin a 

so that 

The maximum ordinate increases more rapidly than the sine of the angle of 

departure. 

Incidentally, multiplying (76) by y' , (82) by x' and adding, we find 

/oA\ ^ X \Xs X ) X \Xs X ) 

(84) p > ,, ^, = 

X + y IT 

From (78) and (60) it is seen that 

The maximum ordinate is increased by a constant following wind. 

In the third of (26), all five factors in the second term to the right are positive, 

whence (x' p)' > — x' n and integrating from to to 4, where to < ts, we find 

x' it(,)p{to) < I x' ndt, 

J to 

since p (4) = by (53), and consequently, since p, = 1 , 

Xs - x{to) 
— ts — to — rrr\ — ' 

X (to) 

Comparing (61) and (19), we therefore see that 

The increase in maximum ordinate dv£ to a constant vertical wind from the time to 
onward is greater than the deflection at the summit due to a constant cross wind of 
the same velocity from the same time to onward. 
Using (76) in (63), we find that 

An increase in the ballistic coefficient increases the maximum ordinate. 
Similarly, it is shown that an increase in atmospheric density ( bH > every- 
where on the trajectory) or in retardation {W > ) will decrease the maximum 
ordinate. We shall finally show that 

The increase in maximum ordinate due to an increase in the ballistic coefficient 
satisfies the inequality 

(85) dys > {2ys — x^ tan a ) -rf ' 

which is interesting inasmuch as the right hand member equals the (approxi- 
mate) expression for hy^ given by the Siacci theory.* 

* See for instance Charbonnier, Balistique exterieure rationnelle, vol. 1 (Paris, Doin, 1907), 
p. 349. 
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To prove this, we observe that in (32), where now k = X = 0, we have 
y' > , x' V + y' p > and ju S 1 , whence 

[x'E{x'v + y'p)]' < - 2Ex'y' 

and integrating from ttota, 

x'E{x' V + y' p) > C'2Ex' y' dt = 1 - x' y' - gx^ = x' y' - g{xe - x) , 

the integration in finite terms being readily verified by means of (1). Con- 
sequently 

Eix'p + y'p)>y' -g^^^, 

X 

pE{x' V + y' p)dt > p(^y' - 9^^)dt 

= \ 2y + {x, — x)^. = 2ys — Xs tsLn a , 
L a; Jo 

and the substitution of this in (63) gives (85). 

6. A SPECIAL INTEGRABLE CASE, AND SOME CONJECTURES REGARDING THE 
RANGE VARIATIONS IN THE GENERAL CASE 

The case to be considered is that of an atmosphere of constant density and 
a retardation proportional to the velocity, so that h = 0, G' = 0. Here 
"K, H,v , p may be expressed rationally in terms oix,y ,x' and y' . Considering 
the range variations only, for the sake of brevity, we have X = 1 , the first 
equation (43) becomes ix' = 0, so that p. = cot w by (44), while the first and 
third of (26) become (x' v)' = — x' and {x' p)' = — x' cot co , so that, by (42) , 

V = -. ) p = -. cot CO . 

X X 

From (45), it is now seen that 



E (x' V + y' p) =x' + y'cotw — g cot u -^—, — , 

X 



whence 



I " E{x' V + y' p)dt = a; + 2t/ cot CO + cota)(a;„ — x)^ 
= a;„ ( 1 — cot CO tan a ) . 
Consequently, formulas (47) to (52) become 
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8x^ 

x^ 


- ( \ I tan a \ 5»o 
\ tan wjvo 


(initial velocity), 


dx^ 

X., 


= ( cot CO — tan a) 5a 


(angle of departure), 


Sx^ 


'. . x„ — x{U)' 

= L- to- , \ '■ Wx 


(following wind). 


5a;„ 


" , , x„ — x{to)l . 

= ta— to TJT^^ cot CO-Wy 

X {to) J 


(vertical wind), 


Sa;„ 

x^ 


_ / ^ tan a \ 6C 
V tan CO J (7 


(ballistic coefficient). 



These expressions lead to the following conjectures in the general case, which 
have been verified on a large number of computed trajectories, and which it 
would be very interesting to prove even under assumptions on G ( « ) and H{y ) 
stronger than (64) : 

(a) The range variation due to an increase in initial velocity satisfies the 
inequality 



foo. / :. , tana\ 5»o 
x^ \ tan 03 J Vo 



(with or without the restriction (71) ? ) . 

(b) For ^t < tu, V is positive, so that, by (49), the range increase due to 
a constant following wind is less than what would be obtained by a rigid dis- 
placement of the original trajectory with the velocity of the wind during the 
time interval from to to t^ . 

(c) The range increase due to a constant vertical wind from the time to 
onward is less than cot co times the deflection at the point of fall due to a con- 
stant cross wind with the same velocity and from the same time to onward. 

(d) The range variation due to an increase in the ballistic coefiicient satisfies 
the inequality 

dx^ -^ / 1 _ tan a \ 5C 
x^ \ tan CO / C 

where the expression to the right is the one occurring in the Siacci theory 
(Charbonnier, I.e.). 

Technical Staff, 

Office of the Chief of Okdnance. 



